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It has been recently found that for the near extremal Kerr black holes appearing of Zero Damped 
Modes (accompanied by qusinormal mode branching) signifies about inapplicability of the regime 
of small perturbations and the onset of turbulence. Here we show that this phenomena is not 
limited by Kerr or Kerr-Newman solutions only, but also takes place for rotating dilatonic black 
holes for which we have found Zero Damped Modes both numerically and analytically. We have 
also shown that, contrary to recent claims, there is no instability of a charged massive scalar held in 
the background of the rotating dilatonic black hole under physically adequate boundary conditions. 

Analytic expression for dominant quasinormal frequencies is deduced in the regime of large coupling 
qQ , where q and Q are the held and black hole charges respectively. 
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I. INTRODUCTION 

A remarkable feature of the characteristic (quasinor¬ 
mal) oscillation spectra of black holes has been recently 
observed in [1]: in the regime of high black hole rota¬ 
tion a kind of the modes branching occurs. One of the 
branches consist of the modes which for the Kerr black 
hole are given by the following formula 

~ y - A (i,m)y/e-i (n+ f) (1) 


where e = y/l — a/M vanishes in the limit of the extreme 
rotation. Thus, when approaching extreme rotation, 
quasinormal modes with vanishing damping rate, that 
is, very long and asymptotically infinite lifetime appear. 
These modes were called Zero Damped Modes (ZDMs) 
and the real oscillation frequencies of these modes were 
found in [ 2 ] under assumption of smallness of their damp¬ 
ing rates. 

Existence of such long living modes which were ob¬ 
tained initially within formalism of small perturbations 
of black holes might mean that regime of small pertur¬ 
bation does not describe near extremal black hole ade¬ 
quately anymore. Indeed, as it has been recently shown 
in Q, ZDMs signify a qualitatively new, turbulent state 
of the Kerr black hole. Moreover, if the black hole is 
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described by the Kerr solution and may achieve near ex¬ 
tremal rotations, observational consequences in the grav¬ 
itational wave structure and accreting processes become 
at least theoretically plausible [3]. 

The modes branching as well as accurate calcula¬ 
tion of ZDMs were extended to Einstein-Maxwell system 
(Kerr-Newman solution) Hi. An appealing question 
is whether the phenomena of Zero Damped Modes (and 
thereby of the black hole turbulence) occurs only for Kerr 
black holes may be answered after understanding the rea¬ 
son inducing such a peculiar behavior of the spectrum. 
Therefore, here we shall try a straightforward way and 
will see if the phenomena of ZDMs in the near extremal 
regime is not limited by the Kerr (or Kerr-Newman) so¬ 
lution, and takes place also for Sen black holes [gj. We 
shall show that the both the bifurcation and ZDMs do 
exist for Sen black hole and, unlike the Kerr-Newman 
case HI, the phenomena take place also for not small 
electric charge Q. 

Another motivation for analysis of quasinormal spec¬ 
tra of the charged scalar field in the rotating dilatonic 
background is a recent statement that the charged field 
is unstable in its background [13]. Here we shall show 
that this claim is not related to fields’s (in)stability in 
its usual sense, but is related to the bound states of the 
field. Here we have also calculated fundamental quasi¬ 
normal frequencies and found analytical formula for uj in 
the regime of large qQ. 

The paper is organized as follows in Sec II we give 
the basic formulas for the dilatonic metric for the Sen 
black hole [6] including the separation of variables for the 
scalar field equation in the black hole background. Sec. 
Ill is devoted to numerical analysis of qusinormal modes 
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with the help of the Frobenius method and deduction of 
the analytic expression for QNMs in the regime of large 
qQ. Sec. IV discusses effects of spectrum bifurcation and 
Zero Damped Modes, including deduction of analytical 
formula for ZDMs. Sec. V comments on superradiance 
and stability of a massive charged scalar field. 


II. SCALAR FIELD IN THE ROTATING 
DILATONIC BLACK HOLE BACKGROUND 

The Sen black hole Q (axion-dilaton black hole with 
zero NUT 0) is given by the following line element [ij 

ds 2 = — a sin 2 Odp) 2 — E + d0 2 ^j (2) 

~ Ao sin 6 2 _|_ a 2 )d(p] 2 , 

where 

A r = r 2 — 2(M — b)r + a 2 , 

E = r 2 + 2 br + a 2 cos 2 0. 


where S{0) obeys the following equation 


( d? d 

+ cot 6 


m 


2 2 - 2 /) 
a uo sm 0 


\ dO 2 dO sin 2 6 

+2mauj + A — p 2 a 2 cos 2 6) S(6) = 0, 


( 7 ) 


and A is the separation constant. 

This equation can be solved numerically for any value 
of uo in the same way as the equation for a massive scalar 
field in the Kerr black-hole background 0- Let us note 
that, when p = 0, Eq. 0 can be reduced to the well- 
known equation for the spheroidal functions. In this case, 
for any fixed value of uj the separation constant A can be 
found numerically using the continued fraction method 
m. When the effective mass is not zero, the separa¬ 
tion constant A(cj,/i) can be expressed, in terms of the 
eigenvalue for spheroidal functions A (cj) 0 , as 


A(cj, p) = A(\/cj 2 — p 2 , 0) + 2 ma{\/ uj 2 — \i 2 — uS) + p 2 a 2 . 


When a = 0, one has A = £(£ +1), £ — 0,1, 2_ For 

nonzero values of a, the separation constant can be enu¬ 
merated by the integer multipole number £ > \m\. 


a is the rotation parameter, M is the ADM mass and 
M > b > 0. The Maxwell and dilaton fields are given by 


III. QUASINORMAL FREQUENCIES AND THE 
FROBENIUS METHOD 


A^dx ^ 1 


Q^=(dt — a sin 2 Odcp ), 

t* 2 + a 2 cos 2 0 

VF---, 


where the electric charge Q is related to b as 


( 3 ) 

( 4 ) 


The radial function R(r) satisfies the following wave 
equation |l3] 

d d K 2 \ 

-A r - + --^ + 2 br) -x) R { r )=0, (8) 


Q 2 = 2 WMb. 


where 


We shall parameterize the metric by the following three 
nonnegative parameters: the event horizon r+, the inner 
horizon r_, and the parameter 6, 

0 < r_ < r + . 

Then, the black hole’s mass and rotation parameter can 
be written as 


K = l jj{r 2 + 2 br + a 2 ) — ma — qQr. 

Equation (j8j) has an irregular singularity at spatial in¬ 
finity and two regular singularities at r = r±. We impose 
the quasinormal boundary conditions, implying purely 
outgoing wave at spatial infinity and purely ingoing wave 
at the event horizon r + . Then, the appropriate Frobenius 
series is 


2 M = r + + r_ + 26, a 2 = r + r_. 

A massive charged scalar field satisfies the 
Gordon equation, 


Dag^V^gDp + n 2 ) V> = 0 


\V~9 


R( r ) — 

Klein- 

where 

K 

(3) a dA r /dr 


( r-r+ 
\r - r_ 


r=r_|_ 


— 2 icy. 

e ip (r — r_) 2cr 7/(r), (9) 


dp 

dr 


K 

A r 


+ — /i 2 CJ, 


where 

d 

Da ^ ~dx^ ~ 

q and p are the field’s charge and mass respectively. 

One can separate variables by the ansatz [9] 

^ = e- iujt + im( t > S(0) R (r ), (6) 


and a can be find from the regularity condition of the 
equation for y at spatial infinity. 

For the QNMs the function y(r) is regular at the hori¬ 
zon and spatial infinity, so that the series in the vicinity 
of the event horizon 

rfr) = (^£+)‘ 

k=0 V 7 
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converge at r = oo 0 . 

The coefficients ak(w) satisfy the three-terms recur¬ 
rence relation, which has the form 

0*0(11 + fiodo = 0 5 

OLkOLk+l + PkCLk + lkdk-l = 0 , k = 1 , 2 ,.... (10) 

The coefficients a^, /3/e, and can be found in a closed 
form. For the massless case (fi = 0) the corresponding 
expressions take the following form 


a k = (k + 1) ^2i(cjr + (r + + r_ + 26) — qQr+ — am) 

— (k + l)(r+ (11) 

/ 3 /e = Aa 2 qQuo - 8r + (6 + r + )(r + + r_ + 26)cj 2 
+z(3r + — r_)(l + 2 k)(qQ — 2buj) 

+8(am + gQr + )(6 + r + )u; (12) 

+2i(l + 2k) (am — 2r^_cj) 

+ (2& 2 + 2& + 1 + A)(r+ — r_) 

+4gQr + (r + + 26 )cj - 4 qQ(am + gQr + ), 


7fc = 


k — 2i(uj(r + + r_ + 26) — gQ) I I fc(r + — r_) 


—2i(u;r + (r + + r_ + 2b) — qQr+ — am) . (13) 


In order to find QNMs we solve numerically the follow¬ 
ing equation (jjl, fig 


0 = /3 0 - 


«o 7 i ^172 <^273 
/3i — /32 — /33— 


(14) 


For the contour plots we calculate logarithm of the ab¬ 
solute value of the righthand side of the above equation 
. For more details on Frobenius method see Sec. 3 of 

In order to find the overtone number n we rewrite <0 
as 



O^n— l7n O^n— 27n—1 O^n— 37n—2 ^07l 

/3n—1 fin—2 fin—3 /3o 

C^n7n+1 ^n+l7n-|-2 ^n+27n+3 

fin-\-l fin-\-2 /3 n +3 — 


(15) 


and search for the most stable root. 

Following 4], we observe that for, large gQ, <a n = 
O(gQ), fin = O(qQ) 2 , and y n = O(gQ) 2 , implying 
that m reads 


Pn 

(qQ ) 2 


= O 



(16) 


Considering A <C ( qQ ) 2 , from ([T6]) we find that 
cj 1 


Vn : lim — =-—. (17) 

qQ^oo qQ r _)_ T r_ T 26 


In this regime, fi n a n = 0(1), fin^> In — 0(1), and, 
as shown in |4j, 


A = 2auo 



i — m 

2 


+ o( 1), 


(18) 


where the brackets denote the integer part. Then (fl5|) is 
reduced to 




(19) 


from which we find the subdominant contributions to the 
QNMs for qQ 1, 


2Muj — c<j(r_|_ T r — T 26) — qQ T — 

r+ 



2n + 1 
4 


TO 



1 

2 


( 20 ) 


The above formula is in excellent agreement with numer¬ 
ical data on fig. [lj 

We have also calculated the fundamental quasinormal 
modes (i = m = 0) for various parameters of rotation as 
a function of 6 (see fig. El. There it can be seen that in¬ 
creasing of the dilaton parameter 6 from 0 up to value of 
a few r + leads to considerable decreasing of the damping 
rate by quite a few times, so that dilatonic black hole has 
much longer living modes for moderate values of 6. Then, 
for larger 6 the damping rate slowly increases, apparently 
achieving some constant asymptotically. The real oscil¬ 
lation frequencies monotonically decrease and also seem 
to approach constant values in the regime of large 6. 


IV. BIFURCATION AND ZERO DAMPED 
MODE: NUMERICAL AND ANALYTICAL 
TREATMENTS 

The smaller the multipole number i the higher is ” criti¬ 
cal” rotation at which the bifurcation takes place. There¬ 
fore, in order to demonstrate the effect, it is reasonable 
to choose some high values of £ and non-zero m. On 
figs. [3j 0]it is shown the logarithm of the absolute value 
of difference between the left and right hand sides of the 
continued fraction equation da as a function of real and 
imaginary parts of u. On figs. [3j 0] the left branches cor¬ 
respond to the so-called Zero-Damped modes (ZDMs) 
which have been found for Kerr black holes in [1]. ZDM- 
branch asymptotically approaches pure real modes in the 
limit of extreme rotation. 

For Kerr-Newman black holes Q, for small values of 
black hole’s charge Q the mode branching occurs owing 
to extremal rotation and not owing to extremal charge. 
Larger values of charge Q break down the bifurcation 
effect. Remarkably, for the Sen black hole, the mode 
branching occurs even for not small values of charge in 
the near extremal limit. 
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FIG. 1: Fundamental quasinormal modes i — mn — 0 for various values of b and a (from top to bottom): b = 0.25r+, a = 0.95r+ 
(blue), b = r+, a — 0.95r + (black), b = 0.5r + , a = 0.75r + (green), b as 2r + , a = 0.5r + (magenta). In the units of mass the 
asymptotical behaviour does not depend on b: blue and black lines approach the same limit. 


Im(6t/)M 



ReMM 



b /r+ 


FIG. 2: Fundamental quasinormal modes £ — m — 0 for various values of a (for large b from top to bottom): a = 0.5r+ (blue), 
a = 0.75r+ (green), a = 0.95r+ (magenta). 


In HI it was shown that ZDMs satisfy a-^0asa^r + . 
Acting in similar fashion with [l|, we conclude that ZDMs 
have the following form 


m±qQ n 

2(r + +b) + 



( 21 ) 


Quasinormal modes depicted on our plots (figs. [3j 0]) is 
in complete concordance with the above expression for uo. 
Non-ZDMs do not scale as oc (r+ + 6) -1 leading to larger 
separation between ZDMs and non-ZDMs for higher val¬ 
ues of b (see Fig. [3j). 


V. REMARK ON STABILITY AND 
SUPERRADIANCE 

Recently, it was claimed in [l3] that a massive charged 
scalar field in the background of Sen black hole has unsta¬ 
ble modes, which might be also important for the black 
hole stability. As previous publications on quasinormal 


modes of charged scalar fields in various asymptotically 
flat or de Sitter black hole backgrounds, including Kerr- 
Newman and dilaton ones showed no unstable modes 
d, IH, [TH, this motivated us to consider here also the 
dominant quasinormal modes in the range of parame¬ 
ters, according to 0, where the instability could be 
expected, namely, the range determined by the super¬ 
radiant regime. 

To show the range of frequencies that gives rise to the 
superradiance effect, we need to study the asymptotic 
behavior of radial solution R(r) at r —> r + and r oc. 
Taking r —>> oo the general solution to ((8]) is a superposi¬ 
tion of ingoing and outgoing waves 

R = exp ir\/uj 2 — fi 2 + 0(ln(r))^ (22) 

+7^exp uj 2 — fi 2 + 0(ln(r ))^ , 

where 1Z is the reflection coefficient. 

There is only ingoing wave near the horizon 

R = (r — r+)~ ia (T +■ 0(r — r+)), (23) 
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FIG. 3: Level plots of the logarithm of the absolute value of the difference between the left and right sides of the equation with 
continued fraction as a function of the real (horizontal axis) and imaginary (vertical axis) parts of the frequency for the scalar 
field i — 10, m = 7, a = 0.985957r+ (r_ = 0.9721117*+), from left to right: b = 0, b — 0.1 r+, b = 0.5r+. 
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FIG. 4: Level plots of the logarithm of the absolute value of the difference between the left and right sides of the equation with 
continued fraction as a function of the real (horizontal axis) and imaginary (vertical axis) parts of the frequency for the scalar 
field i — 10, m = 7, b = r+, from left to right: a = 0.97r+, a m 0.98r+, a = 0.985957r+. 


where T is the transmission coefficient. 

It is possible to prove for real frequencies that TZ > 1 

only if — < 0, in what follows that [13] 


0 < cj < 


qQr+ 


2 Mr, 


or 




2 Mr, 


< cj < 0. (24) 


It is worthwhile of notice that the superradiance condi¬ 
tion for the Sen black hole coincides with the Kerr black 
hole, but differs for the Kerr-Newman black hole (cf., 0)- 
We have made thorough numerical search of quasinor¬ 
mal frequencies in the above range of parameters and 
found no unstable quasinormal modes (an example of 
non-existence of the unstable mode, expected in |13j can 
be seen in fig. [5]). Therefore, no any true instability oc¬ 
curs for the quasinormal spectrum of a massive charged 
scalar field in the Sen background, contrary to [13]. The 
reason of this discrepancy is straightforward: boundary 
conditions used in [l^ is not the one, corresponding to 


black hole’s response to the initial perturbation, that is, 
not the quasinormal mode boundary condition. The lat¬ 
ter implies purely ingoing waves at the event horizon and 
purely outgoing waves at infinity, while the boundary 
conditions in [13[ correspond to bound states of the field 
0 - Instability of such bound states is known for a long 
time j20| , occurs for a massive field in the rotating black 
hole’s background and simply means that the correspond¬ 
ing particle will change its quantum state in order to 
move from the superradiant sector. Thus, bound states 
of a neutral massive field are unstable around Kerr black 
hole in the regime of superradiance [20], but no insta¬ 
bility takes place under quasinormal boundary condition 
0 . Generally speaking, superradiance do can lead to 
the instability in the classical sense, but these are cases 
of asymptotically non-flat backgrounds, such as, for ex¬ 
ample, asymptotically de Sitter background 0- Notice, 
that in the latter case, no rotation is necessary to trigger 
the instability. 
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FIG. 5: i — m = 1, — 0.1, a — b — 0.5r+, gQ — 0.1, 

cjr+ = 0.44148 — 0.07654z is the fundamental mode, u = /jl 
is a local minimum which can be approached as close as we 
want. The requirement uj < /i along the real axis does not 
allow to impose the quasinormal boundary conditions. 

VI. CONCLUSION 

We have shown that the effect of bifurcation of the 
quasinormal spectra and appearance of Zero Damped 


Modes takes place not only for Kerr, but also for Sen 
black holes. Unlike for the Kerr-Newman black hole, even 
highly charged black holes have Zero Damped Modes in 
the limit of extreme rotation. 

This means that inapplicability of the regime of linear 
perturbations for the near extreme black hole’s rotation 
is not a peculiar property of the Kerr solution, but oc¬ 
curs also for other models of rotating axisymmetric black 
holes. Apparently, the effects of ZDMs and the bifur¬ 
cation occurs for much larger class of geometries than 
Kerr-Newman or Sen black holes, and further research 
may find possible connections between these phenomena 
and black hole geometries. 

We have also shown that, contrary to the recent claim 
0, quasinormal spectrum of a charged scalar field in 
the background of the rotating dilatonic black hole does 
not have unstable modes. Analytical expressions for Zero 
Damped Modes (in the near extremal regime) and for 
“ordinary” quasinormal modes (in the regime of high qQ ) 
were deduced. Numerical data for QNMs are in excellent 
agreement with the analytic formulas. 
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